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Abstract 

Nucleation, commonly associated with discontinuous transformations between 
metastable and stable phases, is crucial in fields as diverse as atmospheric science and 
nanoscale electronics. Traditionally, it is considered a microscopic process (at most 
nano- meter), implying the formation of a microscopic nucleus of the stable phase. 
Here we show for the first time, that considering long-range interactions mediated 
by elastic distortions, nucleation can be a macroscopic process, with the size of the 
critical nucleus proportional to the total system size. This provides a new concept 
of "macroscopic barrier-crossing nucleation". We demonstrate the effect in molecular 
dynamics simulations of a model spin-crossover system with two molecular states of 
different sizes, causing elastic distortions. 
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FIG. 1: Nucleation and domain formation for short-range interaction systems, a, 

A schematic example of a droplet of a short-range magnetic interaction system. The blue circle 
with radius r is a cluster of down-spin molecules in an up-spin phase (red part), b, Microscopic 
competition of the surface and bulk free energies of a droplet. The surface free energy is an 
increasing function of the radius of the droplet (r) (upper broken line) and the bulk free energy is 
a decreasing function of r (lower broken line). The solid line is the sum of these energies (AE). 
The critical nucleus with the radius r c gives the maximum of the total free energy, c, Snapshots of 
nucleation for a short-range interaction model (Ising model). Red and blue denote up and down 
spins, respectively 

Nucleation is a barrier-crossing process, 1 3 in which a metastable phase decays via a 
critical nucleus for which the increase in surface free energy is compensated by the bulk 
energy decrease .^E3 if the cluster becomes bigger than the critical size, it grows, while if 
smaller, it shrinks. The size of a critical nucleus is determined by microscopic competition 
between the surface and bulk free energies of a microscopic cluster, and thus the size of the 
critical nucleus is microscopic 3 7 9 (see Fig[l[iandb). To be precise, this situation is realized 
in short-range interaction systems, where separation of the energy between the bulk and 
surface is allowed. Nucleation theoried^^ have been based on this idea, and so far only 
microscopic nucleation is known. 

Consider a typical short-range interaction system at low temperatures: the rf-dimensional 
Ising model 11 defined by the Hamiltonian, % = JJ2i,j — h^i where a = ±1 (up and 
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down spins). The free-energy barrier for a droplet with radius r is AE = —C^h^ + C^Jr^ 1 . 
Here C\>r d is proportional to the volume of the droplet and C s r d ~ l is proportional to the 
area of the phase boundary. As depicted in Fig. [TJ), the radius of the critical droplet is 
given when the droplet has the maximum excess free energy as r c = Cs ^~^ J . The critical 
radius r c is independent of the system size. Figure shows an example of nucleation in a 
circular system (open boundary conditions, OBC) for the Ising model (d — 2). Nucleation 
takes place both in the bulk (inside) and at the boundary. Nucleation at the boundary is 
energetically more favorable, but when the system becomes larger (the bulk-to-boundary 
ratio becomes large), nucleation in the bulk becomes dominant (Supplementary 1). 

However, when the interaction is of long range, the nature of nucleation is different. Be- 
cause it has been pointed out that an elastic interaction due to lattice distortion causes 
an effective long-range interaction,^ the nucleation process in systems with elastic in- 
teractions (e.g., spin-crossover systems j^"^' martensitic system^^Q] anc [ J a hn- Teller sys- 
temd^EHHl) g ] 10U | ( j ^ e investigated. In this work we present properties of the nucleation in 
long-range elastic interaction systems with OBC. 

In molecular crystals, e.g., transition-metal complexes, a molecule often displays bistabil- 
ity in both its electronic state and molecular size (structure). External stimuli, e.g., change 
of temperature, pressure, photoirradiation, etc. change the molecular size. The distortion 
caused by the change of size induces an elastic interaction, which acts as an effective long- 
range interaction. 1 12125 1 26 ^ Spin-crossover (SC) compounds are a typical example of the above 
situation (see Figj2]a-c), where the low-spin (LS) and high-spin (HS) states are separated 
by an energy barrier, and the LS molecule is smaller than the HS one. Indeed, SC systems 
show a wide variety of phase transitions under external stimuli jEHUD 

In the present work we study the nucleation dynamics of circular (d = 2) crystals of 
a long-range elastic interaction system and show that the nucleation is a barrier-crossing 
process. However, the size of the critical nucleus (r c ) is proportional to the system size (R). 
Thus, a macroscopic nucleation mechanism is realized, which is qualitatively different from 
previously known nucleation mechanisms. 

Results 

We adopt the following Hamiltonian for the model 

N 2 N N p2 

*.= gf^g^M + E^ (i) 
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(iJ) 

The first and second terms describe the motion of the intramolecular mode of the zth 
molecule. The radius of the molecule is r^, and the conjugate momentum is p^. The mass 
for the motion is m (Supplymentary 2). The intramolecular potential energy l/ mtra (x^) is 
shown by the solid curve in Fig. [2^i, where Xi = T{ — r^s- Here r^s is the ideal radius of 
the LS molecule, and that of the HS molecule is rns = ^ls + Ar. V^ mtra (^) provides a 
symmetric vibration mode and induces changes of the molecular size. The third and fourth 
terms describe the ceneter-of-mass motion of the molecules [X^ Pi) with mass M. The 
intermolecular interaction V^ ter (see Methods) for nearest and next-nearest neighbors is a 
function not only of the coordinates Xi and Xj, but also of the molecular radii Ti and rj. 
Although the interaction terms in this Hamiltonian appear only to be of short range, elastic 
interactions mediate the effect of the local lattice distortions over long distances.^ We study 
relaxation and nucleation in this model at low temperatures by using a molecular dynamics 
method.^ The critical temperature of the first-order transition of the model is T c ~ 0.9, 
below (above) which the LS (HS) state is the equilibrium stated 

We observe relaxation from the metastable HS state to the LS state at a low temperature 
(T = 0.2) in approximately circular crystals on a square lattice. This temperature is much 
lower than the critical temperature T c . For the initial states of the relaxation (the metastable 
HS phase), we gave a set of velocities to all molecules according to the Maxwell-Boltzmann 
distribution by using a random number sequence. 

Snapshots of the configuration during the course of a relaxation event are depicted in 
Fig. [2]l and Fig. [2^, where the diameters of the circular crystals are 200 and 100 particles 
(denoted as 2R = 200 and 2R = 100), respectively. Figures [2]ll and [2^1 show configurations 
when the HS fraction (/ H sP reaches the value / HS ^ 0.95 for 2R = 200 and 2R = 100, 
respectively. Nucleation starts from one point along the circumference. The subsequent con- 
figurations are given in Figs[2]l(2— 6) and [2^(2— 6) for 2R = 200 and 2R = 100, respectively. 
The corresponding values of /hs are the same in both systems. 

As we show below, the configurations of Figs[2]ll and [2^1 are those of the critical nucleus, 
and Figs [2]l(2— 6) and Figs [2^(2—6) correspond to deterministic growth of the LS droplet 
after the formation of the critical nucleus. It should be noted that the LS domain shapes 
can be well characterized by using the contact angle (wetting angle) 29 of n/2. We checked 
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FIG. 2: Elastic interaction system and nucleation features, a, Intramolecular potential 
energy V(x) shown by the solid (blue) curve. The dotted curves are low spin (LS) and high spin 
(HS) potential energies without quantum mixing. The curvature for the LS state is 4 times larger 
than for the HS stat^l i n this work but other choices for the ratio of the curvatures, for example 
equal curvatures, do not change the essence of the results. The energy unit is 100—300 K for SC 
compound:^ (Supplementary 2), and it is also the unit of the temperature T. b, LS molecule 
(blue) and HS molecule (red). The HS molecule is larger in size, c, Distortion due to the difference 
of the molecular sizes, d, e, Snapshots of the configuration during relaxation from the HS phase 
for (d) 2R = 200 and (e) 2R = 100. The value of / H s is 0.95 for dl and el, 0.90 for d2 and e2, 
0.71 for d3 and e3, 0.50 for d4 and e4, 0.26 for ,d5 and e5, and 0.15 for d6 and e6. 



the configurations for other relaxations from different initial conditions (different random 
number sequences for the molecular velocities) and found the same features of nucleation 
and growth. Here the size of the critical nucleus is found to be proportional to the system 
size. The shapes of the critical nuclei and also the following clusters are geometrically similar 
in systems of different sizes. This feature is qualitatively different from that of short-range 
interaction systems, 30 in which the critical droplet has a specific size independently of the 
system size. 

In Figs [3^i and [3)3 the time dependence of /hs is shown for systems with 2R =100 and 
2R = 200, respectively. Because a single nucleation event dominates the process, the escape 
time from the metastable state is random and governed by a Poisson process. However, once 
nucleation starts, the process is almost deterministic. The crossover between the stochastic 
and deterministic regimes determines the critical nucleus size and a threshold value of /hs- 
These are typical characteristics of barrier-crossing dynamics. 

To capture this feature, we study the relaxation of /hs after passing a given value of 
/hs(= /hs)- The passing time t tT is defined as /^s = /Hs(*tr)- As mentioned above, the 
time evolutions after passing the threshold value, i.e., /hs(* ~~ ^ tr ) are expected to overlap. 
We plotted the data of /hs(* — ^ tr ) f° r various trial values of /jj Sj and found /Hs(^ tr ) — 0.95 
gives the threshold as depicted in Fig. ^jp. The same value is observed in both systems with 
2R = 100 (Fig. [3)3) and 2R = 200 (Fig. [3]i), and we conclude that it is independent of the 
system size. 

Discussion 

To examine the features of the critical nucleus and check the size dependence, we analyze 
the total potential energy of the system (E tot = V^ ntra + V^ nter ) as a function of the relative 
size of the LS domain. It is considered that the entropy effect is small enough compared to 
the energy barrier during the relaxation at this low T. As a parameter to characterize the 
domain size, we define 9 as the central angle. With the contact angle of 7r/2, the domain 
region (lens-shaped part) is defined for any 9 (Fig. |4^l), where the interface between the 
two phases is given by the circle of the radius = i?tan(0/2), whose center is the crossing 
point of the two tangential lines. The value of E tot for a given 9 is obtained as follows. In 
the circle of the HS phase, we replace HS molecules in the lens-shaped part subtended by 
9 by LS molecules. Then we move all molecules slowly so as to reach the minimum total 
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FIG. 3: Relaxation processes from the metastable HS phase, a, b, HS fraction (fas) 
versus time at T = 0.2 for (a) 2R = 100 and (b) 2R = 200. The time of the collapse of the 
metastable state is widely distributed, c, d, The three panels of (c) show the time dependence of 
/hs 5 starting from the time t tr , for 2R — 100 when /^ s =0.98, 0.95, and 0.9 (from left to right). 
The three panels of (d) correspond to the case of 2R = 200. 

potential-energy state, and obtain the energy value of this stationary state. We define the 
energy density as p = ^jp, where TV is the number of molecules in the system and TV ~ ttR 2 , 
and also the relative energy density: Ap = p — po =0 as the difference between p of the 
stationary state and that of the complete HS phase (p$=o). 

We show Ap as a function of 9 for several system sizes (2/2) in Fig. [4^l. For small values of 
9(< 7r/10), Ap is almost constant and then Ap increases with 0. In this region the cluster is 



expected to shrink in the relaxation process. Around 9 = 2.37r/10, Ap shows the maximum 
value and it decreases for larger 9. It should be noted that at this #(=2.37r/10 ) /hs is equal 
to 0.95, which agrees with the threshold value of /hs i n the analysis of the relaxation curves 
(Fig§3 and d). 

When the droplet size exceeds the critical size, the domain expands. For different system 
sizes, this critical size of the droplet (Yd) changes in proportion to the system syze (R). 
Namely, the critical angle exists, but not a specific critical size. This fact was demonstrated 
in Fig. [2]l and Fig. [2^, i.e., the domain shape is almost the same for systems of different size. 
Thus we call this process "macroscopic nucleation" , and we believe that it should hold even 
in the bulk (continuum) limit. We depict Ap at 9 = 2.37r/10 (peak position) as a function 
of l/R in Fig. [ij), and find the dependence: Ap = Ap - The value Ap « 0.035 is 

considered the value of the bulk limit. 

The behavior of "macroscopic nucleation" is qualitatively different from that observed in 
short-range interaction systems. The bulk and surface contributions to the potential barrier 
cannot be distinguished in this long-range interaction system, which is similar to interface 
energies of binary alloys due to elasticity 31 , and the elastic interactions suppress both bulk 
nucleation and multi-droplet nucleation at the boundary (Supplementary 1 and 3). 

In summary, we propose a new concept of "macroscopic nucleation" for systems with 
long-range interactions. The domain formation exhibits geometric similarity for circular 
crystals of any size. This means that the size of the critical nucleus is proportional to the 
system size and macroscopic nucleation is realized. Recognition of this mechanism should 
give important insights for all systems in which local structural changes cause a distortion of 
the lattice. In addition to the spin-crossover type systems considered here, the mechanism 
should hold for martensitic and Jahn- Teller systems, etc. 

Methods 

The role of the intermolecular potential is to release the local distortions due to the difference 
of the sizes of neighboring molecules. For this purpose, we adopt the following potential,^ 
l/^ nter ( X X) Xj,r h rj) = f(dij - Ar), where f(u) = D (>'(") + e -b'(u-u )y The varia bl e 
Uq is a constant such that f(u) has its minimum at u — and = | Xi — Xj \ — (r^ + rj). 
For nearest neighbors, Ar = 0, o! = 0.5 and b' — 1.0 are set, and the energy minimum is 
realized when the neighbors have the same size. For next-nearest neighbors, Ar = 2{y/2 — l)f 
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FIG. 4: Barrier-crossing of macroscopic nucleation. a, The excess energy density Ap as 
a function of 9 for the system sizes 2R = 100, 200, and 300. The inset is the definition of 9 and 
domain region is given for 9 using the contact angle of tv/2. For all R, the values of /hs = 0.98, 
0.97, 0.95, 0.92, 0.90, 0.86, and 0.83 are given at 9/tv = 1/10, 1.7/10, 2.3/10, 2.9/10, 3.5/10, 
4.2/10, and 4.8/10 respectively. After 9/tt = 4.8/10, Ap decreases monotonically until 9/tt = 2 
(LS phase), b, The dependence of Ap on 1/R at 9 = 2.3tt/10 for 2R = 100, 200, 300, 400, and 
500. Ap approaches a finite value as R approaches infinity. 

with f = (tls + ?"hs)/2, a' = 0.1 and b' = 0.2 are set. This provides a small force sufficient 
to ensure the stability of the crystal structure (this is specific to coordination z — 4). 
The parameter D associated with the strength of the intermolecular interaction was set to 
D = 20, which is strong enough to cause a first-order phase transition. Here the type of the 
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potential function is not so important, and the basic mechanism of macroscopic nucleation is 
universal for other types of intermolecular potentials (harmonic or anharmonic potentials). 
The other parameters were set as rns = 9, = 1, and m — M — 1 (Supplementary 2). 
Molecular dynamics simulations were performed using a Nose-Hoover thermostat.^ With 
this method, the timescale of the simulation is influenced by the thermostat parameters. 
Here, we used this effect to our advantage to perform the simulations for large systems in a 
computationally feasible time. Although we sacrifice the ability to measure nucleation times, 
which we expect to increase dramatically with system size, our method allows us to observe 
the scale-invariant spatial structure of the nucleation process in systems of very different 
sizes, as shown in Figs [2]l and [2^. The qualitative feature of nucleation in a stochastic 
Poisson process, followed by deterministic growth shown in Fig. [3| is also preserved. 
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